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Abstract. Nodes in a complex networked system often engage in more than one 
type of interactions among them; they form a multiplex network with multiple 
types of links. In real-world complex systems, a node's degree for one type of 
links and that for the other are not randomly distributed but correlated, which we 
term correlated multiplexity. In this paper we study a simple model of multiplex 
random networks and demonstrate that the correlated multiplexity can drastically 
affect the properties of giant component in the network. Specifically, when the 
degrees of a node for different intera<;tions in a duplex Erdos-Renyi network are 
maximally correlated, the network contains the giant component for any nonzero 
link densities. On the contrary, when the degrees of a node are maximally anti- 
correlated, the emergence of giant component is significantly delayed, yet the 
entire network becomes connected into a single component at a finite link density. 
We also discuss the mixing patterns and the cases with imperfect correlated 
multiplexity. 
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Correlated multiplexity and connectivity of multiplex random networks 



2 



1. Introduction 

In the last decade, network has proved to be a useful framework to model structural 
complexity of complex systems [3 [5] . By abstracting a complex system into nodes 
(constituents) and links (interactions between them), the resulting graph could be 
efficiently treated analytically and numerically, through which a large body of new 
physics of complex systems has been acquired O El [5] . Most studies until recently 
have focused on the properties of isolated, single networks where nodes interact with 
a single type of links. In most, if not all, real-world complex systems, however, nodes 
in the system can engage in more than one type of interactions or links; People 
in a society interact via their friendship, family relationship, and/or more formal 
work-related links, etc. Countries in the global economic system interact via various 
financial and political channels ranging from commodity trade to political alliance. 
Even proteins in a cell participate in multiple layers of interactions and regulations, 
from transcriptional regulations and metabolic synthesis to signaling. Therefore, a 
more complete description of complex systems would be the multiplex network [5] 
with more than one types of links connecting nodes in the network. Multiplexity can 
also have impact on network dynamics; many dynamic processes occurring in complex 
network systems such as behavioral cascade in social networks [7] or dynamics of 
systemic risk in the global economic system f5| should be properly understood from 
the perspective of multiplex network dynamics. Since its introduction in the social 
network literature [5], however, only a handful of earlier related studies have existed 
in the physics literature, notably Refs. OITnilTT], and the understanding of generic 
effects of multiplexity remained lacking. 

More recently, related concepts such as interacting networks 12J and 
interdependent networks fl3j have been introduced and studied. Leicht and D'Souza 
|12| studied what they called interacting networks, in which two networks are coupled 
via inter-network edges, and developed a generating function formalism to study their 
percolation properties. Buldyrev et al. [13 studied the interdependent networks, in 
which mutual connectivity in two network layers plays an important role, and found 
that catastrophic cascades of failure can occur due to the interdependency. Although 
the specific contexts are different in these studies, if one regards each type of links 
in a multiplex network to constitute a network layer, and the multiplex network as 
multilayer network, the multiplex networks and the interacting or interdependent 
networks may be described by a similar framework at the mathematical level. In this 
sense, these studies have provided a pioneering insight relevant to multiplex networks 
that there can be nontrivial effects of having more than one type of links, or channel 
of interactions, in networks [UKTS]. 

In Refs. |12L I13j . network layers were coupled in an uncorrelated way, in the 
sense that the connections or pairings between nodes in different layers are taken 
to be random. In real-world complex systems, however, nonrandom structure in 
network multiplexity can be significant. For example, a person with many links in 
the friendship layer is likely to also have many links in another social network layer, 
being a friendly person. Such a nonrandom, or correlated multiplexity has recently 
been observed in the large-scale social network analysis of online game community 
[16] , in the world trade system [17] , and in transportation network systems [HJ [19] , 
and its impact on network robustness has been studied |18l I20 | . In this paper, our 
main goal is to understand the generic role of correlated multiplexity in multiplex 
system's connectivity structure. 
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Figure 1. Schematic illustration of multiplex ER networks with three kinds of 
multiplex couplings discussed in the text. MP (MN) stands for maximally-positive 
(maximally-negative) correlated multiplexity. 

2. Model and formalism 

To study how the network connectivity is affected by correlated multiplexity, we 
consider the following model of multiplex networks with two layers, or rfiipZex networks. 
The network has N nodes connected by two kinds of links, modeling, for example, 
individuals participating in two different social interaction channels. We refer the 
subnetwork formed by each kind of links to as the network layer. Each network layer 
I {I = 1,2) is specified by the intralayer degree distribution 7r^'^(fc;), where degree fc; 
is the number of links within the specific layer I of a node. The complete multiplex 
network can be specified by the joint distribution n(fci, fc2) or the conditional degree 
distribution 11(^2 Generalization into I > 2 layers is straightforward. 

The total degree of a node in the multiplex network is given by k — ki + k2 — ko, 
where ko denotes the number of overlapped links in the two layers, which can be 
neglected in the N ^ oo limit for random, sparse networks with largest degree 
of at most order 0{^fN) [3T]. One can obtain the total degree distribution P(k) 
from the joint degree distribution or the conditional degree distribution as P{k) = 
Sfci fea n(fci, fc2)^fc,fci+fc2 ~ X^fci n(fc — fci|fci)7r(^)(fci), where 5 denotes Kronecker delta 
symbol. From P(fc), one can follow the standard generating function technique [22| to 
study the network structure: The generating function gg i^) of the degree distribution 
P{k) of the mutiplex network can be written as 

oo 

go{x) = ^(^)^' = E n(fci, ^2)0;'=^+'=^ . (1) 

A;— ki.k2 

Emergence of the giant component spanning a finite fraction of the network signals 
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the establishment of connectivity. Size of the giant component S is obtained via 

oo 

5 = l-goH = l-^P(fc)ti', (2) 

fc=0 

where u is the smallest root of the equation x = g'Q{x) / §'^{1) = gi{x), that is, 

oo 

Mean size of the component to which a randomly chosen vertex belongs, (s) , plays 
the role of susceptibility and is given by 

is) = l + ^^ (4) 

^ go{u)[l - g[{u)] 

The condition for existence of the giant component (that is, 5 > 0) is given by the 
existence of a nontrivial solution it < 1, leading to the so-called Molloy-Reed criterion 



k 



k{k ~ 2)P{k) ^ {k"^) - 2{k) > . (5) 



It is worthwhile to note the related recent generalizations of the generating function 
method for interacting and interdependent networks [23'. 

For a multiplex network system, the total degree distribution P(fc) is determined 
from the joint degree distribution Il{ki, ^2), which depends on the pattern of correlated 
multiplexity. Therefore, the presence of correlated multiplexity affects the multiplex 
system's connectivity (figure 1). In the following, we specifically consider duplex 
networks of two Erdos-Renyi (ER) layers [26] and three limiting cases of correlated 
multiplexity. Using analytical treatment with mean-field-like approximation as well as 
extensive numerical simulations, we present how the correlated multiplexity can affect 
the emergence of the giant component in the multiplex system. Similar procedures 
can also be applied to multiplex scale-free network models |27) . 



3. Degree distributions 

3.1. Uncorrelated multiplexity 

In the absence of correlation between network layers, the joint degree distribution 
factorizes, n„„corr(fci, ^2) = tt^^-' (fci)7r'^^' (^2). The total degree distribution of 
the multiplex network is given by the convolution of 7r*-'-'(fc/), Puncorr{k) = 
^^^^Q7r(^)(/ci)7r(^)(fc — fci), and its generating function (7q"™'''"(x) = gQ^\x)gj^\x), 

where go \x) is the generating function of 7r('^(fc;). Using 50(2^) — e^'^^^^^ for the ER 
network with mean degree z, we have g^^'^"'^'^ (x) = e(^i+^^)*^^~^) for the duplex ER 
network with mean intralayer degrees z\ and 02, which is nothing but the generating 
function of an ER network with mean degree z\ + Z2. Therefore, we have 



P{k) = (6) 



with z = zi + Z2. 
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Figure 2. (a) Schematic illustrations of how the mean-field- like calculation for 
degree distribution is done, (b, c) Degree distributions of interlaced multiplex ER 
networks with 21 = 22 = 0.7 (b) and 21 = 22 = 1.4 (c). Points are numerical 
simulation results for the uncorrelated (o), MP (□), and MN (o) cases, together 
with lines representing predictions of mean-field-like calculations. 



3.2. Maximally-positive correlated multiplexity 

In the maximally-positive (MP) correlated multiplex node's degrees in different 

layers are maximally correlated in their degree order; the node that is hub in one layer 
is also the hub in the other layers, and the node that has the smallest degree in one layer 
also has the smallest degree in other layers. To obtain the total degree distribution 
P{k) of the multiplex network, we use the following mean-field-like scheme ignoring 
fluctuations in the thermodynamic limit {N — > 00), illustrated for the duplex case in 
figure 2(a). We partition the unit interval into bins of sizes Tr^^'>{ke) sorted in order of 
ke for each £. Combining the two partitions, we have a new partition which can be 
used to reconstruct P{k) as illustrated in figure 2(a). 

3.3. Maximally-negative correlated multiplexity 

In the maximally-negative (MN) correlated multiplex case, a node's degrees in different 
layers are maximally anti-correlated in their degree order; a node that is hub in one 
layer has the smallest degree in the other layer, etc. The mean-field-like scheme for 
duplex case proceeds in a similar way as the MP case, except that we use two partitions 
sorted in opposite orders respectively as illustrated in figure 2(a). 
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Figure 3. (a) The giant component size S and (b) the susceptibility (s) as a 
function of the single-layer mean degree zi of the two-layer multiplex ER networks 
with equal mean layer degrees for uncorrelated (black) , MP (blue) , and MN (red) 
cases. Lines represent the solutions of (2) and (4) under the mean-field scheme 
and symbols represent the numerical simulation results obtained from networks of 
size N = 10* with 10* different configurations for the uncorrelated (o), MP (□), 
and MN (o) cases. Errorbars denote standard deviations. Gray shade denotes the 
region in which 5 = 1 for the MN case (zi > z* = 1.14619322...). 



4. Duplex ER networks with equal link densities 

In this section we consider ER networks with two layers of equal link densities, for 
which the degree distributions are most easily calculated. 

4-1. Uncorrelated case 

The uncorrelated multiplex ER network is simply another ER network with mean 
degree z — 2zi, where zi is the mean degree of the first network layer. The joint 
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degree distribution factorizes, thus the conditional degree distribution n(/c2|fci) is 
independent of ki and n(A:2|fci) = 7r(fc2), where 7r(/c) denotes the Poisson distribution 
with mean degree zi . As we increase zi , we have a second-order percolation transition 
at the critical mean intralayer degree Zc — 1/2 [221 US], at which the giant component 
emerges in the duplex network. Giant component size scales in the vicinity of Zc as 
S {zi — Zc)^ with /3 = 1, and the susceptibility as (s) ~ \zi — Zc\~'^ with 7=1, 
following the standard mean- field critical behaviors |22) . 



4.2. MP case 

In this case, in the mean-ficld-likc scheme {N 00) each node has exactly the same 
degrees in both layers, so the conditional degree distribution becomes n(fc2|A:i) = 
Sk2,ki- Thus we have the degree distribution of the duplex network as 

P(k)-{ e-^^%'^V(fc/2)! (fceven), , . 

^^^^ - 1 (fc odd). 

Therefore, the Molloy-Reed criterion is fulfilled for all nonzero Zi, as (fc^) — 2(fc) = 
4(zi -I- zf) — 2(2zi) = izf > for zi 7^ 0. This means that the giant component exists 
for any nonzero link density, that is, 

= 0. (8) 

In fact, one can obtain the solution of S and (s) explicitly in this case: As P{k) = 
for odd k, (3) has u = as a nontrivial solution, from which it follows from (2) 

S ^1- P{0) = l-e-'\ (9) 

which grows linearly with the link density near origin as S ^ zi, that is /? = 1 . From 
u = 0, the susceptibility (s) — 1 for all zi > 0. This means that only isolated nodes are 
outside the giant component and all the linked nodes form a single giant component. 
All these predictions are confirmed by numerical simulations (figure 3). 



4.3. MN case 

In this case one can easily show that distinct regimes appear as Zi increases. Among 
them, three regimes are of relevance for the giant component properties: i) For 
< zi < ln2, more than half of nodes are of degree zero in each layer so every 
linked node in one layer is coupled with a degree-0 node in the other layer under MN 
coupling. In this regime the conditional degree distribution takes a rather complicated 
form 

r [2^(0) - 1] A(0) (fc2 = 0,fci=0), 
n(fc2|fci) = <^ ^(A:2)/7r(0) (fc2 ^ 0, fci = 0), (10) 

[ 4„o ih ^ 0), 

and thus P{k) is given by 



27r(fc) (fc > 1). 

In this regime there is no giant component, ii ) For In 2 < zi < z* , similar consideration 
leads to n(fc2|fci) and P{k) given by 

(fc2 = 0), 

n(fc2|0)-<( [27r(0) + 7r(l) - 1] /7r(0) (fca = 1), (12a) 
7r(fc2)A(0) (fe2 > 2), 
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Figure 4. Data collapse of (a) scaled giant component sizes, SNP/'', and (b) 
scaled susceptibility, {s)/N^/'', vs. (2 - z^''^)N^/'' for the MN case, with /3 = 1 
and u = 3, consistent with the conventional mean-field behaivor /3 = 1 and 7 = 1. 
N is 10*{A), 10''{n), and 10'^{o). 



[27r(0) + 7r(l) - 1] /7r(l) (^2 = 0), 
n(fc2|l)-<( [1 - 2^(0)] /7r(l) (fc2 = l), (126) 

(fc2 > 2), 

n(fc2|fci) = 4.,o (fci>2), (12c) 

and 

r (fc-o), 

P(u^-) 2[2^(0)+^(l)-l] (fc = l), 

^^''>-\ 2^(2) - 27r(0) + 1 (fc = 2), ^^"^^ 
[ 27r(/c) (fc > 3). 

In this regime, (fc^) — 2(fc) = 2(z^ — zi — 2e~^i + 1), which becomes positive for 
zi > where 

^ = 0.838587497... (14) 

Therefore the giant component emerges at a significantly higher link density than the 
uncorrelated multiplex case. Being delayed in its birth, the giant component grows 
more abruptly once formed [see figure 3(a)] than the other two cases. This regime is 
terminated at z = z* , determined by the condition 27r(0) + 7r(l) = 1, from which we 
have z* = 1.14619322... Hi) For zi > z*, we have F(0) = P(l) = 0. In that case, 
we have u — from (3) and thereby S — 1 from (2). This means that the entire 
network becomes connected into a single component at this finite link density, which 
can never be achieved for ordinary ER networks. Despite these abnormal behaviors 
and apparent differences in the steepness of the order parameter curve near Zc, the 
critical behavior is found to be consistent with that of standard mean-field, /3 — 1 and 
7 = 1 (figure 4). 

5. Imperfect correlated multiplexity 



In the previous section we have seen that maximally correlated or anti-correlated 
multiplexity affects the onset of emergence of giant component in multiplex ER 
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Figure 5. Plot of equation (15) for the critical mean degree Zc as a function of 
g, the fraction of correlated multipex nodes. The cases q = 1 denote maximally 
correlated multiplexity, < 5 < 1 partially correlated multiplexity, and g = 
uncorrelated multiplextiy. The vertical dotted line is drawn atg = 2 — l/ln2 
across which Zc takes different formulae in equation (15b). 



networks. Despite its mathematical simplicity and tractability, in real- world multiplex 
systems the correlated multiplexity would hardly be maximal. Therefore it is 
informative to see how the results obtained for the maximally correlated multiplexity 
are interpolated when the system possesses partially correlated multiplexity. 

To this end, we consider duplex ER networks where a fraction q of nodes are 
maximally correlated multiplex while the rest fraction 1 — q are randomly multiplex. 
Then the degree distribution of the interlaced network is given by Ppartiai{k) = 
QPmaximai{k) + (1 ^ q) Puncorr{k) , wherc maximal is either MP or MN. Following 
similar steps as in the previous section we obtain the critical link density as a function 
of q as 

z,^{l-q)/2 (15a) 

for positively correlated case and 

1/(2 -g) (9<2-l/ln2), 

zi{q) (9 > 2-1/ In 2) ^'^'^ 

for negatively correlated case, where Zi{q) is the solution of {2 — q)zi— zi — 2qe~^^ +q = 
0. In figure 5, we show the plot of Zc as a function of q given by (15). This result shows 
that the effect of correlated multiplexity is not only present for maximally correlated 
cases but for general q. 

6. Duplex ER networks with general link densities 



In the previous sections we focused on the cases with zi = Z2- In this section 
we consider general duplex ER networks with z\ 7^ z^- We performed numerical 
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Figure 6. Numerical simulation results of the size of giant component of duplex 
ER networks with (a) MP, (b) uncorrelated, and (c) MN cases. In (d-f), the 
giant component size S (red) is plotted for 22 = 0.4, along with the assortativity 
coefficient r (blue) for the MP (d), uncorrelated (e), and MN (f) cases. Mean- 
field-like calculation results (lines) deviates from the numerical simulation results 
(o) when assortativity coefficient (□) becomes nonzero, that is, the network 
display degree-degree correlation. Errorbars denote standard deviations from 10* 
independent runs. 



simulations with N = 10^ for MP, uncorrelated, and MN correlated multiplex cases 
with Zi and Z2 in the range from to 3. In figure 6, the giant component size S for 
general zi and Z2 are shown. Similarly to the zi = zi cases, the giant component 
emerges at lower link densities for the MP case but grows more slowly than the 
uncorrelated case, whereas it emerges at higher link densities for the MN case but 
grows more abruptly and connects all the nodes in the network at finite link density. 
Therefore the effect of correlated multiplexity is qualitatively the same and generic. 

It is worthwhile to note, however, that although the mean-field-like approximation 
for P(fc) introduced in section 3, which was very effective for equal link densities, still 
yields a qualitatively correct picture, it fails in quantitative agreement with numerical 
simulation results for general z\ ^ zi- To understand the origin of this discrepancy, we 
consider network correlations. Assortativity coefficient r [28] defined as the Pearson 
correlation coefficient between the (total) degrees of nodes connected by a link in the 
network measures the degree-degree correlations in the network at the two-node level. 
Nonzero assortativity is known to alter the connectivity properties of networks |28) . 
To see the role of degree-degree correlations induced by correlated multiplexity, in 
figure 6 (d-f) we plot the giant component size obtained from numerical simulations 
and mean-field calculations, together with the assortativity coefficient as a function zi 
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with Z2 — 0.4 fixed. With uncorrelated multiplexity, the degree-degree correlation is 
absent, r = 0, and the numerical simulation and mean-field calculation agree perfectly 
[figure 6(e)]. For MP and MN cases, however, the assortativity of multiplex network is 
generically nonzero (positive), that is the multiplex network becomes correlated, which 
is responsible for the deviations between the numerical simulation and mean-field 
calculation results. This discrepancy vanishes at zi = Z2, at which the assortativity 
also vanishes [figures 6(d,f)]. Note that there was no degree-degree correlations at the 
individual network layer level. Therefore this result shows that generically correlated 
multiplexity not only modulate P{k) but also introduce higher-order correlations in 
the multiplex network structure. Meanwhile, it is worthwhile to note that such a 
multiplexity-induced degree correlation has similar origin as the correlation in colored- 
edge networks recently studied in a different context of network clustering [29] . 

7. Conclusion 

In conclusion, we have studied the effect of correlated multiplexity on the structural 
properties of multiplex network system, a better representation of most real-world 
complex systems than the single, or simplex, network. We have demonstrated that 
the correlated multiplexity can dramatically change the giant component properties. 
With positively correlated multiplexity, the giant component emerges at a much lower 
critical link density, which even approaches to zero for MP case, than for uncorrelated 
multiplex cases. Once formed, however, the giant component grows much more 
gradually. With negatively correlated multiplexity, the giant component emerges at 
a much higher critical density than for uncorrelated multiplex cases, but once formed 
it grows more abruptly and can establish the full connectivity to connect the entire 
network into a single component at a finite link density. These results show that a 
multiplex complex system can exhibit structural properties that cannot be represented 
by its individual network layer's properties alone, the impact of which on network 
dynamics is to be explored in future study. 
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